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The  singular  perturbation  method  for  the  supercavitating  propeller  had  four  scaling 
parameters;  1)  span  length  R,  2)  chord  length  c,  3)  blade  spacing  d  and  4)  cavity 
length  The  first  problem  solved  here  assumed  that  c/R  and  ?C/R  were  of  order 

but  d/R  was  of  order  of  unity.  The  nature  of  the  singular  perturbation  problem  for 
such  a  case  was  similar  to  that  for  the  subcavitating  propeller  solved  by  Brockett 
except  for  the  solution  of  the  inner  region.  The  thrust  and  torque  coefficients 
were  obtained  explicitly  without  solving  the  integral  equations.  Since  the  non¬ 
linear  supercavitating  flow  theory  was  employed  in  the  present  work  as  the  inner 
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solution,  there  existed  no  limitation  for  the  flow  incidence  angles  or  blade  pro¬ 
file  shapes  and  thus  the  present  solution  would  provide  more  accurate  results  than 
those  with  the  linearized  theory. 

The  second  problem  treated  here  was  the  case  in  which  c/R  was  of  order  of  e  but 
ic/R  and  d/R  were  of  order  of  unity.  This  was  the  case  having  long  cavities  behind 
the  propeller  blades  so  that  even  when  the  chord  shrank  to  a  line,  the  cavities  were 
left  behind  the  lifting  lines.  This  portion  of  cavity  sneets  was  called  "source 
sheets",  the  singularity  strengths  of  which  were  obtained  through  the  cavity  sheet 
matching, a  totally  different  matching  procedure  from  the  regular  matching.  The 
first-order  inner  solution  used  a  closure  condition,  i.e.,  the  total  source  term  Sq 
equal  to  zero.  As  the  result  of  the  cavity  sheet  matching,  however,  it  was  dis¬ 
covered  that  the  second-order  inner  solution  had  to  use  a  finite  Sq,  the  value  of 
which  was  determined  through  the  matching.  It  seemed  that  the  second-order  matching 
automatically  corrected  the  overstretching  assumption  made  in  the  first-order  inner 
solution. 

Two  other  problems  were  posed,  both  having  small  blade  spacings.  The  outer  solu¬ 
tions  may  be  quite  different  from  those  in  the  above  two  cases  since  if  one  looks 
at  such  a  propeller  from  the  far  field,  the  blade  elements  will  not  be  identified. 
This  may  require  the  actuator  disc  concept  for  the  outer  solution,  with  a  con¬ 
ventional  pressure  jump  across  the  disc  if  the  cavity  is  short  but  with  a  cavity 
pressure  drop  to  be  applied  if  the  cavity  is  long,  similar  to  the  theory  of  Tulin 
(1965) .  These  two  problems  have  not  been  carried  out  here  due  to  the  enormous 
amount  of  work  required  even  for  the  first  two  problems. 
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1.0 


BACKGROUND 


Solving  the  propeller  flow  problem  is  inherently  difficult  due 
to  its  geometric  complexity  as  well  as  the  propeller  rotational 
motion.  For  fully  wetted  propellers,  many  different  types  of 
theories  have  been  developed  to  date,  both  for  design  and  off- 
design  analyses.  These  include  lifting  line  and  lifting  sur¬ 
face  theories,  using  the  singularity  distribution  method.  In 
the  performance  prediction  method  the  strengths  of  sources  and 
vortices  represent  the  thickness  and  loading  of  blade,  respec¬ 
tively  and  are  to  be  determined  by  satisfying  the  f low-tangency 
condition  on  the  propeller  blade  surface.  The  numerical  analysis 
using  high  speed  computer  is  essential  and  sometimes  encounters 
the  instability  problem. 

The  problem  of  supercavitating  (s/c)  propeller  flow  is  even 
more  difficult  due  to  the  additional  feature,  i.e.,  the  existence 
of  the  cavity  on  the  suction  side  of  the  blade.  The  extent  of 
cavity  varies  significantly,  depending  upon  the  ship  speed  and/ 
or  propeller  rotational  speed.  At  the  design  point,  the  super¬ 
cavitating  propeller  usually  operates  with  relatively  short  to 
medium  length  cavities.  However,  as  the  flow  speed  or  rotational 
speed  increases,  the  cavity  becomes  longer  with  increase  in 
thickness.  The  flow  passage  between  blades  is  gradually  blocked 
and  finally  totally  choked  over  the  entire  blade  span.  Although 
the  effect  of  the  thick  cavity  may  theoretically  be  represented 
by  the  source  distribution,  a  difficult  problem  exists  in  that 
the  location  of  the  cavity  is  not  known  until  the  complete 
supercavitating  propeller  flow  is  solved. 

Three  different  types  of  concepts  exist  in  accounting  for  this 
thick  and  long  cavity  effect  on  the  s/c  propeller  flow.  When 
the  sectional  loading  used  for  propeller  analysis  is  calculated 
based  on  the  single  foil  configuration,  the  'retarding'  flow 
correction  due  to  the  cavity  blocking  by  Tulin  (1965)  may  be 
necessary  in  addition  to  the  conventional  induced  flow  correction. 
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This  concept  is  similar  to  that  of  the  actuator  disk  theory, 
except  that  the  cavity  pressure  instead  of  fully  recovered 
pressure  at  downstream  infinity  is  used  in  the  momentum  theory. 
The  effective  incoming  flow  velocity,  different  from  that  of 
the  upstream  infinity,  is  used  as  part  of  the  correction  to  the 
propeller  flow  diagram.  This  method  was  actually  applied  to 
the  recent  Hydronautics  design  of  a  supercavitating  propeller, 
designated  as  Model  7607.02  (see  the  report  by  Bohn  and  Altman 
(1976)).  Experimental  results  for  the  propeller  obtained  by 
Bohn  (1977)  and  Peck  (1977)  showed  that  the  design  method  over¬ 
predicted  the  thrust  by  about  10  percent. 

Instead  of  correcting  the  effective  incoming  flow  with  the 
retarding  flow  effect,  the  second  category  of  the  concept 
utilizes  the  forces  calculated  from  the  cascade  flow  analysis 
for  the  propeller  sectional  loading.  This  type  of  method  is 
widely  used  in  the  design  and  analysis  for  the  pump  and  turbine 
in  which  the  flow  is  usually  well  confined  in  the  circular 
duct.  The  force  coefficients  obtained  in  this  theory  are 
generally  much  lower  than  those  of  single  foil  cases  because 
the  pressure  on  the  blade  is  substantially  lowered  by  the 
cavity  attached  to  the  adjacent  blade.  These  forces  are  now 
used  for  the  induced  flow  corrections  due  to  the  free  vortex 
sheet.  It  must  be  noted  that  the  retarding  flow  correction  is 
not  necessary  because  it  is  considered  that  the  cavity  blockage 
effect  is  already  taken  care  of  when  obtaining  the  sectional 
forces.  The  author  recently  developed  such  a  theory  based  on 
the  above  concept  (see  Furuya  (1976,  1978)) ,  The  prediction 
capability  of  the  theory  is  good  for  low  advance  speed  range 
J  where  a  strong  cavity  choking  condition  prevails.  Inciden¬ 
tally,  an  accurate  prediction  of  the  s/c  propeller  performance 
at  low  J's  is  important  since  the  high  speed  taking-off  condi¬ 
tion  occurs  under  such  conditions.  However,  for  larger  J's  for 
which  the  cavity  has  a  shorter  length,  the  theory  substantially 
deviates  from  the  experimental  data.  It  means  that  the  super¬ 
cavitating  propeller  theory  based  on  the  cascade  data  will  be 
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suitable  for  predicting  the  off-design  performance  at  the  same 
as  or  lower  J  than  of  design  point  but  that  the  accuracy 
for  larger  J's  is  questionable. 

One  of  the  major  questions  arising  from  applying  the  two- 
dimensional  cascade  flow  to  the  three-dimensional  propeller 
flow  is  as  to  what  flow  incidence  angle  must  be  used  for 
determining  the  force.  Since  in  the  cascade  flow,  unlike  the 
single  foil  flow  case,  the  flow  field  is  completely  separated 
by  the  cascade  blade,  there  exist  two  reference  flow  angles, 
one  at  upstream  infinity  and  the  other  at  downstream  infinity. 
In  the  pump  and  turbine  flow  analysis  it  is  conventional  to 
use  the  geometric  mean  flow  angle.  This  may  be  understandable 
since  in  such  devices  the  upstream  flow  field  is  physically 
separated  from  the  downstream  flow  field.  In  the  propeller 
flow  the  flow  field  is  so  called  "singly"  connected  like  the 
single  foil  case.  The  only  difference  is  that  the  cascade 
effect  locally  exists  in  the  former.  It  is,  therefore,  not 
quite  clear  which  angle  is  to  be  used  for  the  present  propel¬ 
ler  analysis,  the  upstream  flow  angle  or  geometric  mean  flow 
angle  of  cascade  theory.  The  above  propeller  theory  of  Furuva 
(1978)  used  the  former  for  the  two-dimensional  cascade  flow 
analysis.  It  is  considered  that  an  appropriate  flow  angle  to 
be  used  may  be  somewhere  in  between  the  above  two  extreme  cases 

The  third  category  is  that  of  Yim  (1978)  who  developed  a  design 
method  by  combining  the  lifting  surface  theory  with  cavity 
effects.  The  strength  of  source  as  well  as  that  of  doublet 
is  determined  by  satisfying  the  boundary  conditions  both  on 
the  solid  and  cavity  boundaries.  However,  the  method  was  devel 
oped  as  a  design  tool,  not  applicable  to  the  off-design  per¬ 
formance  prediction  of  supercavitating  propellers. 

In  the  situation  just  described,  it  is  considered  useful  to 
investigate  the  problem  from  a  different  point  of  view,  namely, 
with  the  singular  perturbation  method  (SPM) .  Brockett  (1972) 
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first  applied  the  method  to  the  subcavitating  propeller  flow. 

He  used  the  linearized  propeller  theory  for  the  "outer"  solu¬ 
tion  and  a  two-dimensional  linearized  theory  of  the  fully 
wetted  single  foil  flow  for  the  "inner"  solution.  Through  the 
matching  procedure  the  circulation  distribution  of  the  outer 
solution  was  explicitly  determined  with  the  induced  flow  cor¬ 
rections  included.  However,  the  scaling  parameters  used  in  the 
work  of  Brockett  were  limited  to  two,  i.e.,  chord  and  span 
lengths,  just  like  the  finite  span  wing  case,  so  that  the  blade 
aspect  ratio  £  was  only  the  small  parameter  for  expansions  to 
be  made.  It  means  that  the  results  of  Brockett  can  only  be 
applied  to  such  propellers  having  large  aspect  ratio  with  small 
number  of  blades,  e.g.,  two  to  three. 

The  singular  perturbation  problem  for  the  supercavi tating  pro¬ 
peller  is  much  more  complex  in  that  there  exist  four  scaling 
parameters.  These  include  1)  span  length,  2)  chord  length, 

3)  blade  spacing  and  4)  cavity  length.  Although  the  blade 
spacing  must  have  been  included  even  in  the  SPM  of  Brockett 
(1972)  for  the  conventional  propeller,  it  was  just  omitted 
there  under  the  condition  that  the  propeller  blades  be  sparsely 
distributed  as  mentioned  above.  Different  combinations  of 
the  above  four  parameters  will  provide  various  types  of  prob¬ 
lems  for  the  supercavitating  propeller  with  SPM.  In  the 
present  study  four  most  typical  cases  are  identified,  as  will 
be  seen  in  Section  3  and  the  solution  methods  for  the  first 
two  problems  among  others  will  be  presented  herein. 

We  will  start  with  the  basic  theories  which  will  form  the 
basis  for  the  inner  and  outer  solutions  (Section  2) ,  which 
will  be  followed  by  the  classification  of  the  problem  (Section 
3)  and  then  its  solution  methods. 
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2.0 


BASIC  THEORIES 


The  singular  perturbation  method  for  supercavitating  propellers 
will  employ  various  two-  and  three-dimensional  flow  theories  as 
inner  and  outer  solutions.  Most  of  these  theories  have  been 
well  established  and  wil^  be  found  much  in  literature.  One 
should  study,  however,  on  which  theory  is  suitable  for  the 
inner  or  outer  solution  and  for  the  matching  procedure.  As  a 
preparatory  work  for  the  major  part  of  the  present  study,  the 
basic  two-  and  three-dimensional  theories  to  be  used  herein  will 
be  described  in  the  following. 
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2 . 1 


LINEARIZED  SUPERCAVITATING  PROPELLER  THEORY 


There  exist  various  ways  of  deriving  the  potential  function  ; 
for  propeller  flows  but  only  a  limited  number  of  papers  are 
available  for  the  supercavitating  propeller  flow.  The  paper 
bv  Cox  (1968)  represents  such  a  theory  and  this  method  will 
be  employed  throughout  the  present  work  as  the  outer  solution. 

Figure  2.1.1  shows  a  schematic  diagram  of  a  supercavitating  pro¬ 
peller  blade  rotating  around  the  x-axis  in  the  uniform  flow  of 
velocity  U  approaching  from  the  negative  infinity  of  x.  The 
center  of  the  propeller  shaft  is  placed  along  the  x-axis. 
Although  the  cylindrical  coordinate  system  (x,  r,  9)  will  be 
used  in  the  present  propeller  analysis,  it  will  be  transferred, 
if  necessary,  into  the  cartesian  coordinate  system  shown  in 
Figure  2.1.1.  The  velocity  components  in  the  x,  r,  9  directions 
are  first  written  in  terms  of  the  perturbation  quantities  for 
U  as  (U  +  u  ,  u  .  u„  )  where  the  subscript  0  is  used  for  the 

X  T  -  Q 

inertial  coordinate  9  in  order  to  distinguish  it  from  the 

rotating  coordinate.  Defining  u  as  (u  ,  u  ,  ua  ) ,  the 

x  r  “  q 

linearized  equation  of  motion  can  be  written 


3t 


vV±  = 

3X 


which  will  be  expressed  in  terms  of  the  velocity  potential  <p 
as  follows 


2 i 

it 


+ 


-P 


(2.1.1) 


where 

u  =  7p  (2.1.2) 

P  =  P/p- 

Solution  for  the  partial  differential  equation  (2.1.1)  is  obtained 
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X  -  X 


(x ,  r ,  3 q  ;  t)  = 


-  if  p(v  +  «-.r.Vt-  (2.1.3) 


where  the  boundary  condition  of  p  =  0  at  x  =  — »  has  been  used. 
Applying  the  continuity  equation  7y]£=0  to  Equations  (2.1.1) 
and  (2.1.2),  the  Laplace  equation  is  obtained  for  p 


'2  o  =  0 . 


12.1.4) 


Use  of  the  Green's  identity  for  p  will  provide 


o  (x 


(sb  0 


dS 


-lit 


Ey.+  S. 


_ _ l 

R„.  3n 

ro 


0  L 


p(x',r',9';t) 


dS 


12.1.5) 


where 

n 

R.„ 


normal  to  linearized  surface,  positive  from  pressure 
side  to  suction  side  . 


(x-x1)2  +  r,2  +  r2  -2r'r  cos  '•?, 


-0 


"  ’o  _  “o 


+  'J 


=  9  '  -  ut 
=  ¥  (k-i} 


K 

S 


§b'§< 


=  number  of  blades 

=  moving  surface  consisting  of  propeller  blade  sur¬ 
face  and  cavity  surface 

=  propeller  blade  and  cavity  surfaces,  respectively. 


The  moving  blade  surface  S  is  considered  to  be  composed  of 
helical  lines  having  varying  pitch  in  the  r-direction;  i.e., 

S:  x-«RMr)  =  0 

•\  (r)  =  U  (r) /wR 


ii 
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where  ,\  denotes  the  advance  coefficient  as  a  function  of  r. 
The  directional  derivative  for  this  surface  is  then  given 


r  -J-  -  r  (0  +  u)t)  R.\  -i-  -  —  — U 

5  3x  0  r  3r  r  36  Q 

— =  rnr 


3^ 


"  2  2  21  1 
r  +  (RAr  +  ( r 6 Q R a r )  ^  2 


where  \  =  d A/dr .  Thus, 


r  (x  -  x  ' )  +  R\r  sin  -  r  O  1  ^  +  uit)  Ra^  (r  '  -  r  cos  j'  ! 


r  +  (RA)^  +  -|r(6,Q  +  u;t)RA|i 


7R3 
—  f 


and 


dS 


d9  dr . 


For  the  steady  state  problem  presently  investigated,  time  t  is 
set  to  zero  in  the  above  equations.  It  means  that  the  rotating 
axis  and  inertial  axis  coincide  instantaneously,  then 

i  (x,  r,  9)  =  t>v  +  } 


where 


t>v  =  potential  due  to  the  vortex 


■-ik //**'•  Hi  ? 

rh  9L  ~  v 


dv 


dS ’dr ’ 


12.1.6) 


=  potential  due  to  the  source 

.  1  3E  r  x_e'An 

■-iSj I 

v-  r\  J 


rh  eL 


d0  '  dr  ' 


(2.1.7) 


r,  =  radial  coordinate  corresponding  to  the  propeller  hub 
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L'  T'  'E 


angular  coordinates  at  the  blade  leading  edge,  trail¬ 
ing  edge,  and  at  the  cavity  end  point,  respectively 


dp  = 


s  = 


P 

M 


p  -  p  (pressure  doublet) 
(pressure  source) 


-%  +  , 

-IE  _  d£ 

9n  9n 


7  2  2 

r  +  (r’0  ’X  ) 


N  = 


D  = 


R  a 


r  1  v  +  r.\  sin  5 - —  (v  -  x)  v  (r  1  -  r  cos  :) 

r 


- r-  ( v  -  x)  -9  +  J, 

V  x 


2  2  2 

v  +  r  •  +  r  -2r  ’  r  cos  p 


-  1 
2 


v  =  dummy  variable  for  integrals, 


It  must  be  mentioned  also  that  the  normalization  has  been  made 
for  all  quantities  above  as  follows, 


(x,  r,  x',  r',  v)  5  (x,  r,  x',  r',  v  ’ )  /R 


d  =  f/RU 

N  =  N  /R2 

v  y 

R  =  R  /R 

>  v 

“r  =  VR 

s  =  s  /2"U2/R 
P  P 

ip  =  ip/2"TU2 


where  R  is  the  propeller  radius.  After  normalization,  the 
bar  above  each  character  has  been  deleted.  Therefore,  it  should 
be  understood  that  all  mathematical  expressions  in  equations 
hereafter,  including  Equations  (2.1.6)  and  (2.1.7),  are  those 
which  were  normalized. 


The  potential  i  due  to  the  source  effect  in  Equation  (2.1.7)  is 
rather  complex.  A  simpler  form  can  be  obtained  if  the  source 
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distribution  is  used  instead  of  the  oressure  source  S  , 

l  e  ‘  p 

K  r  r E 

*  = _ L  V'  /  /  S  (r  '  ,  e  ' )  d9  '  dr  ' 


k  r  r  & 

vlr  f f 

k=l  y  r  9  ' 


(2.1.8; 


where 


,  =  [(x  -  e  ' \)  2  +  r'2  +  r2  -  2r’r  cos  r.  .]  I 


'V  —  a  i  _  o  .  -T 

'9  '  -  0  +  Jk. 


In  Equation  (2.1.8),  s (r 1 , 6  1 )  is  the  source  term  (not  a  pressure 
source  in  this  na  ^  tajH  i  rh  mnroc  4-  +-Ka  ki  _ _ i  j _ 


2.2  2-D  NONLINEAR  SUPERCAVITATING  FOIL  THEORY 

In  order  to  provide  a  high  accuracy  of  the  results  obtained  from 
the  present  SPM,  a  nonlinear  theory  for  the  supercavitating 
(s/c)  foil  will  be  used  as  an  inner  solution.  Many  flow  models 
having  various  cavity  closure  conditions  are  available.  These 
include  the  open  wake  model  of  Wu  (1962)  ,  and  single  and 
double  spiral  vortex  models  of  Tulin  (1964).  Among  others, 
the  single  spiral  vortex  model  of  Tulin  will  be  used,  the  reason 
for  this  selection  being  explained  later. 

The  single  spiral  vortex  model  was  first  employed  by  Larock  and 
Street  (1965,  1968)  in  the  nonlinear  theory,  but  only  applied 
for  the  calculation  of  the  s/c  flat-plate  foil  characteristics 
and  for  the  inverse  specification  of  foil  profile.  The  theory 
and  computer  program  were  recently  developed  for  general  pro¬ 
file  cases  by  Furuya  and  Maekawa  (1980)  and  compared  with 
experimental  data.  The  key  features  of  the  theory  necessary  for 
the  singular  perturbation  method  will  be  summarized  here.  The 
flow  configuration  and  boundary  conditions  are  shown  in  Figure 
2.2.1(a)  in  which  the  physical  coordinate  system  employs  (x,  y) 
in  tune  with  that  used  in  the  inner  region  as  will  be  seen  later. 
The  physical  flow  field  is  mapped  onto  the  potential  plane  W  as 
shown  in  Figure  2.2.1(b)  which  is  then  transformed  onto  the  upper 
half  of  a  new  plane  ;  =  ^  +  in  of  Figure  2.2.1(c)  by  a  mapping 
function 


or 


W 


*1 


(2.2.1) 


(2.2.2) 


The  cavity  end  point  is  mapped  to  infinity  and  the  infinity 
point  in  ~  =  x  +  iy  or  W  plane  is  now  mapped  onto  a  point  ;  =  ia. 
A  hodograph  variable  w  is  introduced 
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where  q  and  9  are  the  magnitude  and  direction  of  the  flow 
velocity.  Thus 


where  q^  and  qc  denote  the  uniform  flow  velocity  at  the  upstream 
infinity  and  that  on  the  cavity  wall,  respectively. 

The  boundary  conditions  on  the  real-axis  *  are  now  expressed 
either  in  terms  of  9  or  x; 

(i)  x  =  0  ,  -  oo  <  £  <  _i  ancj  b  <£  <», 

(ii)  a  =  3  ,  -1<5<0,  S-tan-MOyOO/dxl 

(iii)  9  =  1  +  3,  0  <  5  < 

where  y(x)  denotes  the  foil  profile  shape  of  the  wetted  portion. 
This  is  a  typical  mixed-type  boundary  value  problem,  and  the 
solution  for  to  is  readily  written 


We  have  a  total  of  four  unknown  quantities,  a,  b,  P,  $  , 

Li 

requiring  four  equations  to  determine  them  uniquely. 


The  various  boundary  conditions  will  be  applied; 


(i) 


at 


infinity 
oj  (ia)  =  i in 


2  equations 


(ii)  length  of  arc  =  S  arc:  1  equation 

(iii)  body-cavity  system  closure  condition  (see  Larock  and 
Street  (1965)), 


The  report  of  Furuya  and  Maekawa  (1980)  describes  the  method 
of  solving  the  above  system  of  nonlinear  integral  equations 
for  the  four  solution  parameters  in  detail. 


It  must  be  pointed  out  here  that  not  all  of  the  nonlinear  theories 
quoted  before  are  applicable  to  the  inner  solution  of  the  problem. 
We  should  exclude  those  theories  whose  body-cavity  systems  do 
not  close  in  the  potential  plane  W.  This  point  will  become  clearer 
if  W(v)  is  expanded  as  7 -■*  .  Let's  use  the  double-spiral  vortex 
model  as  an  example,  the  flow  configuration  of  which  is  shown 
in  Figure  2.2.2(a).  Mapping  the  potential  plane  W  =  J+ii 
onto  the  upper  half  of  the  ;-plane,  the  hodograph  solution  ~ 
can  be  obtained  as 
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1 _ 

i  + '  * )  (b  -  rn 


ln\  1  +  j 


i  41:  _  f  i 

+  1)U  *  -  b)  ’  ’  "  ^  t/IT^+TTTT7  -  b 


(2.2.5) 


where 


qe  =  qxe 


(2.2.6) 


thus 


w  =  9  +  i  2n  yq1 


W  =  Ac;  ;  A  =  solution  parameter,  yet  to  be  known. 


From  Equation  (2.2.6), 


1  iw  dW  , 
dY  =  — —  e  xr  dc, 


1  .iut-lf,  .  :  .  1.  i^"(»)  -  -0’  (»)‘  1  . 

d  f  —  — -  e  1  +  iu)  (°°)t  +  - t -  — ~  +  •  •  •  •  ^  As  d^ 

-»  *  -  * 


as  :-*■». 


(2.2.7) 


Integration  of  both  sides  of  Equation  (2.2.7)  provides 


__  ^ 

Y  -s»AKq  (;  +  2  C  +  2  in  (  +  ■  •  •)  ,  as  £ 


(2.2.8) 


where 


Kq  =  e 


=  iu  1  (•*>) 


K2  =  j  iV  (»)  -  .o'  (~)  2jy4. 

Inverting  Equation  (2.2.7),  one  obtains  j 

W  =  A;2  ~  -2  K1  ~  AK  2  +  ‘  ‘  '  as  2  00  • 

The  first,  second  and  third  terms  in  expansion  represent  the 
uniform  flow,  corner  flow  and  source-circulation  flows, 
respectively.  Each  term,  except  for  the  second  one,  can  pos¬ 
sibly  find  its  counterpart  from  the  outer  solution.  However, 
the  second  item,  i.e.,  corner  flow  part,  cannot  be  matched 
with  the  outer  solution.  This  illustrates  the  difficulty  of 
using  any  flow  model  which  does  not  close  in  the  downstream 
infinity,  particularly  in  the  potential  plane.  Note  the  dif¬ 
ference  of  the  potential  planes  in  Figure  2.2.1(b)  and  Figure 
2.2.2(b);  one  is  closed  and  the  other  is  open.  It  is  for  this 
reason  that  we  have  chosen  a  single  spiral  vortex  model  for 
the  inner  solution. 


1 


¥ 


I 


( 
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2.  3 


2-D  LINEARIZED  SUPERCAVITATING  FLOW  THEORY 


♦ 


* 


We  will  utilize  the  results  of  the  nonlinear  theory  as  the 
inner  solution  for  the  singular  perturbation  method.  As  will 
be  seen  later,  however,  a  linearized  theory  is  necessary  for 
conveniently  determining  the  source  strength  on  the  body- 
cavity  system.  The  most  suitable  linearized  theory  for  deter¬ 
mining  such  source  strength  is  that  of  the  singularity  distri¬ 
bution  method.  Since  the  first  introduction  of  the  method  by 
Tulin  (1953) ,  many  researchers  presented  papers  for  various 
types  of  problems.  However,  the  basic  concept  is  described  in 
the  original  paper  by  Tulin  for  the  drag  body  case  and  in  that 
of  Fahner  and  Spiegel  (1966)  for  the  lifting  body  case.  The 
linearized  flow  configuration  is  shown  in  Figure  2.3.1  and  the 
linearized  boundary  conditions  for  the  perturbed  velocities  u 
and  v  are  given  as  follows: 


0 

“  2 

*i 

f 

0  <  x  <1, 

+ 

0 

11 

(>1 

0 

=  2 

*i 

/ 

1  <  *  <  v 

1 

O 

II 

f  kN 

=  qx- 

f  ’  (x) 

/ 

0  <  x  <1, 

y  =  0, 

where 


If  the  source  s^  (x)  (0  <  x  <  1)  ,  S2  (x)  (1  <  x  <  9,^)  and  dipole  y  (x) 
are  distributed  along  the  slit  shown  in  Figure  2.3.1,  i.e., 
for  0<x  <2.  ,  the  induced  velocities  due  to  these  singularities 
will  be  given 


+ 

u  = 


+ 


u 


u 


u 


(x ' ) 
x  -  y 


dx ' 


0 


_1_ 

2t 


c 


S2  ( x_' ) 

X  -  X  ' 


dx ' 


1 


0  <  x  c  1 


(2.3.1) 


1  <  x  <2. 

c 


(2.3.2) 
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-1 


+ 

V 


Sj  (x) 
2 


V 


S1  (x) 
2 


0  <  x  <  1 


+  _  S2U) 
2 


v 


s  2  ( x ) 
2 


1  <  x  <  L 


0 


0  <  x  <  1 


0  <  x  <  1 . 


(2.3.3) 

(2.3.4) 

(2.3.5) 

(2.3.6) 


The  boundary  conditions  on  the  slit  are  therefore  satisfied  by 
combining  Equations  (2.3.1)  to  (2.3.6),  providing  a  system  of 
integral  equations  with  respect  to  s.(x),  s2 (x)  and  y(x), 


-L  f  fill!  dj.  +  i  f  till.. 

2n  J  x  -  x'  x  2ir  J  x  -  x' 


1  /  s2  ^  —  Y (x)  a  — 

1  ‘  dx'  +  -Ap-  =  j  qx,  0  <  x  <  1 

(2.3.7) 


1  /'si(x')  -  i  r~  s2(x')  _ 

tc  /  — — =7-  dx '  +  ■=—  /  — — =7-  dx'=^-q,  ,1<X<1 

2  J  x  -  x '  2tt  y  x  -  x’  7 


(2.3.8) 


1 

2tt 


/*  Y(x') 

J  x' 


dx  1 2 


s^x) 


q^ •  f '  ( x)  ,  O^x  <1.  (2.3.9) 


Applying  the  inversion  formula  given  in  Appendix  A  to  Equation 
(2.3.8),  the  following  result  will  be  obtained 


(2.3.10) 


where  a  singularity  at  x  =  i  has  been  placed  for  representing 
an  abrupt  termination  of  the  cavity.  s2(x)  is  now  substituted 
into  Equation  (2.3.7), 
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(2.3.11) 


Application  of  the  inversion  formula  again  to  Equation 
yields  an  expression  for  s^(x)  in  terms  of  7 (x)  ; 


(2.3.11) 


(2.3.12) 


Substitution  of  s^(x)  in  Equation  (2.3.12)  into  Equation  (2.3.9) 
eliminates  s^(x),  and  results  in  the  following  integral  equation 
for  y (x)  only; 


(2.3.13) 


Change  of  variable  with 


(2.3.14) 


will  transform  Equation  (2.3.13)  into  the  following  form 


0  <  <  <  e  (2.3.15) 


where 


,2  ,2 

l  <  _  <.  <  ' 

X(<)  =  — - =-  ,  X-  (<  '  )  =  —2 - Y 


1  +  < 


1  + 


(2.3.16) 


e  =  1/  (  ic  -  1) 
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Finally,  the  dipole  distribution  of  y(x)  is  obtained  by  invert¬ 
ing  Equation  (2.3.15); 


(x(k))  =  (cos  3) 


1  1 
,  ,  2  /e  -  <\ 2 


Y  (sin4 


+  < • COS  J 


1  _e 


(■*.■)  PH5/ 


^ ,  0., 

,  .  ,2  <  -  <  ' 


1  +  < 


(2.3.17) 


where 


3  =  tan  ^  (e) 


and  the  leading  edge  singularity  has  been  introduced  above  for 
representing  the  stagnation  flow  within  the  framework  of  linearized 
theory . 

s^(x)  and  s2(x)  will  then  be  calculated  by  substituting  Equation 
(2.3.17)  into  (2.3.12),  followed  by  (2.3.12)  into  (2.3.10) 


!g(x(<))  =  j  qx  (cos  S)  2  ^  2  ^sin1^  ~ 


vU-2)  H* 


i  e  l  _ 

\  2  f  l  <’  \  2  f  '  (  x  '  ( <  ' 

vy  i-w  i+.-2 


'  ))  d<  ' 
2  <  +  < 


~q  x  f '  ( x  ( < )) 


(2.3.18) 


s2(  x(<)) 


j q^ (cos  S ) 


P 

sin  2 


3 

cos  2 
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3.0 


SINGULAR  PERTURBATION  METHOD  FOR  SUPERCAVITATING 
PROPELLERS 


Application  of  the  singular  perturbation  method  to  the  super- 
cavitating  propeller  flow  will  lead  to  multiple  problems.  Most 
of  the  hydrodynamic  problems  solved  with  this  method  have  only 
two  characteristic  lengths  in  the  flow  field.  Three-dimensional 
fully  wetted  and  cavitating  single  foil  flows  were  solved  by 
many  researchers  including  Van  Dyke  (1964),  Shen  and  Ogilvie 
(1972),  Leehey  (1973),  Furuya  (1975)  and  James  (1975).  The 
two  characteristic  scales  for  such  flows  were  chord  and  span 
lengths  with  its  ratio  assumed  to  be  small.  The  problem  of 
thin  foils  with  rounded  leading  edges  was  also  treated  with 
the  singular  perturbation  method,  for  the  fully  wetted  flow 
by  Van  Dyke  (1964)  and  for  the  supercavitating  flow  by  Furuya 
and  Acosta  (1973).  Again,  two  characteristic  lengths,  i.e., 
chord  and  leading  edge  radius,  existed. 

For  the  supercavitating  propeller  flow,  however,  we  have  four 
characteristic  lengths.  These  are 

i)  propeller  radius  (R) , 

ii)  blade  chord  length  (c) , 

iii)  cavity  length  (lc),  and 

iv)  propeller  blade  spacing  (d) . 

With  these  four  physical  parameters,  many  combinations  of 
scaling  can  be  considered  and  these  are  shown  in  Table  3.0.1. 
Schematic  flow  configurations  are  also  given  in  Figure  3.0.1. 

Case  I  presents  a  problem  which  is  exactly  the  same  as  that  of 
fully  wetted  propeller  flows  but  with  a  difference  existing  in 
the  flow  configuration  of  the  inner  solution.  Instead  of  using 
a  fully  wetted  single  foil  flow,  a  supercavitating  flow  must  be 
employed  for  the  inner  solution.  In  this  case,  the  ratios  of  c  to  < 
and  d  to  R  are  assumed  to  be  all  of  the  same  order  £  so  that,  as 
the  propeller  blade  collapses  to  a  line,  all  pertinent  informa- 
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tion  for  the  propeller  flow  can  be  expressed  around  the  lift¬ 
ing  line. 

However,  this  is  not  the  case  for  Case  II  in  Table  3.0.1.  Since 
we  have  a  long  cavity  behind  each  blade,  even  with  the  blade 
shrinking  to  a  lifting  line,  the  cavity  portion  fails  to  col¬ 
lapse  to  a  line  but  remains  as  a  finite  length  in  the  outer 
flow.  This  portion  cannot  possibly  be  represented  by  any 
quantities  on  the  lifting  line  but  requires  the  concept  of 
source  sheet. 

It  will  be  readily  understood  that  the  strength  of  these  source 
distributions  cannot  be  matched  through  the  expansion  of  the 
inner  and  outer  solutions  around  the  lifting  line.  The  reason 
is  that  the  source  strengths  of  such  cavity  sheet  have  nothing 
to  do  with  the  scaling  of  c/R  =  0(e)  but  belong  to  the  other 
scaling,  i.e.,  2  /R.  It  is  considered  that  the  determination 
of  the  source  sheet  characteristics  requires  a  different 
type  of  matching  with  the  inner  solution  as  will  be  seen  later. 

The  type  of  problems  associated  with  Cases  III  and  IV  are  more 
complex.  As  the  chord  shrinks  to  a  line  for  the  outer  solution, 
the  spacings  between  the  blades  will  become  unrecognizable 
since  the  spacing- to-radius  ratio  is  assumed  to  be  of  the  same 
order  as  that  of  the  chord-to- radius  ratio.  The  outer  solution 
will  have  to  be  solved  with  the  actuator  disk  concept.  If  the 
line  source  term  SQ  is  equal  to  zero  and  no  source  sheet  exists 
(most  likely  in  Case  III),  only  a  careful  consideration  for 
choice  of  the  inner  flow  configuration  may  be  necessary.  How¬ 
ever,  if  Sq  is  not  zero  followed  by  a  source  sheet,  the 
actuator  disk  will  have  to  add  a  source  disk  as  well  as  the 
lower  pressure  region  behind  it.  Thi.  ill  account  for  the 
flow  retardation  effect  in  the  heavily  cavitating  propeller 
flow,  first  explained  by  Tulin  (1965). 

In  the  present  work,  due  to  the  limited  time  available,  only 
the  first  two  cases  will  be  solved  with  the  singular  perturba- 


tion  method,  leaving  the  other  two  cases  to  the  near  future 
effort  when  resources  become  available. 


SHORT  CAVITY  PROPELLER  CASE  (CASE  1) 


3.1.1  Outer  Solution 


An  outer  solution,  valid  in  the  outer  region,  will  be  obtained 
as  a  limiting  form  of  p  when  the  blade  shrinks  to  a  lifting 
line.  First  look  at  j>  in  Equation  (2.1.6).  By  defining  the 
third  integral  in  (2.1.6), 


x-0  1  A 


I(r' ,9') 


/N 


(3.1.1) 


we  can  expand  I(r',9')  for  small  9',  i.e.. 


I(r', 3’)  =  I  (r  '  ,  0)  +  9’  +  H.O.T. 


(3.1.2) 


e  '=o 


where  H.O.T.  stands  for  the  neglected  higher-order  terms, 
0  (9 ' 2) ,  etc. ,  and 


I (r ' ,0)  = 


(3.1.3) 


X-0  '  V 


9  I  ! 
3  9  '  1 


3  [  Nv 

597  R,3 


r'x  +  r A  sin  4* , 


dv  =  -  A 


3.1.4) 


2  '2  2 

Rq  =  x  +r  +r  -2rr’  cos  (-6  +  o^) 


(3.1.5) 


\  =  -0  +  °k- 


We  now  have  in  expansion 

v  a 


1  "  T 


-  -I  E 


U! 

rh  eL 


Ap (r  ' , 9  ' ) 


x 

fb 


r  '  x  +  r  A  sin  4\ 


dv  -  A  9 


+  H.O.T.  dr 'd9  '  . 


(3.1.6) 
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Redefining  some  of  the  terms  in  this  integral,  i.e., 


,,  (r.e,i£i±i!)2 


ip(r,e: 


j 

f 


:(r')  =  /  Yc(r,,e,)d0 


“lx(r'>  E 


T 

J  X0  '  Ip (r  '  ,  0 


’ )  r  ’  d0  ' 


■/ 


Uly(r,)  =  /  1^9 ' Ap (r  '  , 9  ' )  d9 '  , 


where  the  normalization  factor  for  F(r),  i.e.,  2tRU,  has  been 
used  and  y  denotes  the  vortex  distribution. 


Equation  (3.1.6)  can  be  written  as  follows 


= 


k=l 


dv 


1 

xy.  (r1 )  +  r  sin(-9  +  a.  )  *  u.  (r'J 
- — - - - - - =2 -  dr  '  +  H.O.T. 

R0 

It  is  seen  from  Equation  (3.1.7)  that  the  potential  $  is  com¬ 
posed  of  a  line  circulation  as  the  lowest  term,  moment  terms 
having  the  axes  in  the  direction  of  chord  and  that  normal  to 
the  blade,  and  H.O.T. 


(3.1.7) 
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Similarly,  the  potential  ?  in  Equation  (2.1.8)  will  be  expanded 
for  small  8'  with  the  present  assumption  of  short  cavity  in 
mind.  As  8  ' ->-0 , 


(3.1.8) 


and 


(3.1.9) 


therefore , 


9')  de  ' 


(-  Ar  *  x  +  r  X  sin  t1^)  d9'+  H.O.T. 

(3.1.10) 


With  a  new  definition  for  the  following  integrals 


Vr'> 


l 

■/ 


s(r\S')  d6'  , 


3lx<r 


slr'.e'le1 ( — \ r ’ )  dS' 


Sly<r'> 


■/ 


s ( r ' , 9 ' ) 9  '  d9 ' 


Equation  (3.1.10)  will  be  written  as 


■  -5  E  [/ 
k=i  L< 


s0(r,) 


dr 


'*/ 


x  s,  (r  '  )  +  r  sin  3,  •  s.  (r  '  ) 
lx  k  ly 

7372 


dr 


+  H.O.T. 


(3.1.11) 


In  this  equation  SQ/  s^  and  s2x  can  be  interpreted  as  a  line 
source,  the  first  moment  of  source  in  the  direction  of  x  and 
y,  respectively. 


For  matching  purposes,  the  asymptotic  limit  of  the  outer  solu- 

2  2  2 

tion  t>  as  sQ  -  0  is  required  where  s^  =  x  +  y  .  The  first 
term  in  Equation  (3.1.7)  represents  the  potential,  denoted  by 
due  to  the  vortex  line  and  its  free  vortices.  It  is 

V  i 

rather  difficult  to  expand  this  term  directly  and  thus  the 

velocity  components,  u  and  u. will  be  used  for  this  pur- 

ai  c. 

pose,  i . e . , 


u  „ 
ai 


3o 


vr 


3x 


u 


tr 


(3.1.12) 


where 


'V 


—  /> 
k=l  J  p  R , 


(3.1.13) 


R 

- »  j 


Taking  derivatives  of  with  respect  to  x  and  9,  we  find 


u  „  =  u  __  +  u  __ 
a:  a. B  a. F 


(3.1.14) 


ut:  ut:B  + 


Ut:F 


where 


u 


aFB 


k=l 


1 


r ( r ' ) dr ' 


r  sin  V, 
k 


(3.1.15) 


(3.1.16) 
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(3.1.17) 

(3.1.18) 


X (r-r  ' 


cos  7  ) 


+  {x  -  X  (x  +  9)}  r  ' 


sin  y 

T 


dIJ 

(3.1.19) 


=  —  P  4-  n 

k  '  +  °K 


’f'  =  r  +  a. 


■  [<> 


Xt)  -  2rr  ’  cos  V  + 


2  , 2 
r  +  r 


1/2 


(3.1.20) 


and  the  subscripts  B  and  F  for  the  velocity  components  denote 
"Bound"  and  "Free"  vortices,  respectively. 


The  dominating  term  in  u  of  Equation  (3.1.16)  comes  from 

a- a  2  2  2 

k=l,  the  vortex  line  of  its  own,  as  =  x  +  y  -*  0 ,  therefore 

1 


u  ^  _  r_sLn±  f 
aTB  2  J 


r(r' )dr' 


+  H.O.T. 


(3.1.21) 


R 


00 


where 


00 


2  ,  ,  ,  2 
.  +  (r-r ' ) 


2  2  2 
s0  =  x  +y 


y 


r  sin  9 . 


The  integral  in  Equation  (3.1.21)  is  singular  as  s^  -  0  and 
cannot  be  expanded  easily  as  it  is.  One  of  the  techniques, 
among  others,  to  resolve  this  type  of  situation  is  an  appli¬ 
cation  of  the  Fourier  integral  representation  for  1/R  ^  ,  the 
method  having  been  used  by  Ogilvie  (1972)  and  James  (1975)  in 
similar  problems,  i.e., 

_3_ 

1  _  1  3  /  1  \  _  /  2  ^  A  2  - 

r030  "  S°  ^so  '  00/  ~  V°  a/ 

00 

=  -  —  A  '  —  /  KA(jm;sA)  exp(-ima)  dm  (3.1.22) 

so  aso  rJ  0  0 

■"  —CO  J 


where  KQ  is  the  modified  Bessel  function  and 


a  =  r  -  r 


(3.1.23) 


The  expansion  of  KQ  (  j  m  |  sQ )  can  be  made  as  sQ  0  (see  Handbook 
of  Mathematical  Functions  by  Abramowitz  (1964)), 


K.  (  jm  I  sj  - 


2 

-  In  sQ+0  (s Q  In  sQ) 


(3.1.24) 


Substitution  of  (3.1.24)  into  (3.1.21)  and  (3.1.22)  yields 


aTB  2 


/:<r')dr'^/-  ^  (-:n 


SQ  +H.0.T.)  • 


exp ( -ima) dm 


i 

r  sin  0  f 


f(r')dr '/  exp(-ima)  dm  +  H.O.T. 


/*[/■■ 


(r')exp  j-im  (r-r')|  -  dr' 
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+  J  .  (r  ' )  exp  j-i  m  (r  1  -r)}  •  dr  ' 


+  H.O.T. 


:  3 .  i .  2  5 ; 


Integration  by  parts  is  applied  to  Equation  (3.1.25), 

r 


u. 


r  sin  6  f  ,  [21  (r)  /*  d  ?  ( r  ’  )  exp  j-im  ( r-r  1  )[ 

fB~  2^s  2  J  mLim  J  dr'  lm 

0  rh 


dr  ' 


Jm 

[  dr<r' 

7  dr' 


)_  exp  j-im  (r  1  -r ) 
-i  m 


dr1 


+  H.O.T. 


(3.1.26) 


Now  the  method  of  contour  integral  will  be  used  by  defining 


/exp  (- 
1 


exp (-ip^  Z 1 


dA  =  0 


(3.1.27) 


where  the  contour  C  to  be  taken  is  shown  in  Figure  3.1.1  and 
p^  =  r-r'>0.  Thus, 

-IT 


/exp (-if 
i  m 


-ip.,  m)  r  exp(-ip.  (Rcos  5  +  Ri  sin  9)  ... 

-  dm  +  t  +  /  - - - -  i  Re  '  d  9  =  0 


iReJ 


Since  the  last  integral  is  zero  as  R  -*■  =°, 

00 

exp (-ip^m) 


j  exp  ( 


dm=  -71 ,  o=r-r'>0 
m  ‘1 


(3. 1. 28  ) 


Similarly 

exp ( —  ip 2  m) 


/exp ( - 


m 


=  +  tt,  p_  =  r  1  -r  >  0 . 


3.1.29) 


u  in  Equation  (3.1.14)  becomes 
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> 


» 


u 


a7  B 


r  r  1 

r  sin  9  f  dT(r')  ,  ,  /*  d7(r' 

*77^  [7  ‘-,ldr  7  — 


)  _ 


dr  1 


+  H.O.T, 


-  r  (r) 


r  sm 


>-  S0 


+  H.O.T, 


-  r  ( r ) 


x2  +  y2 


+  H.O.T. 


(3.1.30) 


7  {  v*  \ 

where  7(1)  =  7(r.)  =  0  and  /  -1-: - dm  =  0  have  been  used. 


m 


In  the  similar  manner,  u  in  Equation  (3.1.18)  can  be  expanded 

+  H.O.T. |  (3.1. 31) 


Vs'71'1 


-  x 


[_x2  +  y2 


It  is  well  known  that  the  induced  velocities  due  to  the  free 

vortex  sheets,  u  and  u.  in  Equations  (3.1.17)  and  (3.1.18),  have 
a  2  2Z 

finite  values  as  x  +  y  -  0,  i.e.,  on  the  lifting  line  itself. 

The  final  forms  of  the  expansions  for  u  _  and  u  _  are  there- 
fore  written 


u  ~  7  ( r )  —5-^ — zr  +  u  +  H.O.T. 

3.1  Z  Z  cl 

x  +  y 


(3.1.32) 


Ut7-’:(r)-TZ2L2 
x  +  y 


+  u  +  H.O.T. 


(3.1.33) 


3y  using  the  following  relationship 


-1  v  y  -x 

tan  =  —z~ - zr  1  +  — = - zr  1 

x  2  2  —  2  2  J- 

x  +  y  x  +  y 


(3.1.34) 


;  „  becomes 

V  i 


. _  —  7  ( r)  tan  1  + 


_vr  - - -  x  ■  (ua* X +  Ufc- y)  +  H.O.T., 


2  2  . 
a  s  x  +  v  -  0  , 


(3.1.35) 
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The  moment  term  of  in  Equation  (3.1.7),  i.e., 

V  ^  VM 

be  expanded  by  using  the  method  of  Ogilvie  (1970) 


VM 


X 


xu^x<r'  )  +  r  sin(-3)  j^tr') 


dr  ' 


can  also 


uix(r),x+  Vr)  ~y 


2  2 
x  +  y 


+  H.O.T. 


(3.1.36) 


where  only  k  =  1  has  been  used  for  the  leading  ;°rm  of  ; 


VM' 


For  the  potential  due  to  the  source,  again  the  method  of  the 
Fourier  representation  for  1/Rq  and  that  of  Ogilvie  are  used 
for  the  expansion,  which  provides 


o  -  Sn  (r)  1 n  (x2  +  y2)  2 


ilx,r,-X+SlV(r)'y 

2  2 
x  +  y 


+  H.O.T. 


3.1.37) 


Combining  ?v„,  4>  and  in  Equations  (  3.1.35  ,  36  and  37), 
the  two-term  expansions  of  the  outer  solution  <»>  can  finally 
be  written  as 

1 

3  -  Jq  +  7  (r)  tan  1  +  SQ  (r)  in  (x2  +  y2)2  +  (u^  (r)  x  +  ufc  (r)  •  y) 


ulx(r)  -x  +  ijy  ( r )  -y  s^r)  -x  +  sly(r)  -y 


2  ,  2 
x  +  y 


2  2 
x  +  y 


H.O.T.  (3.  1. 38) 


where  5n  is  the  uniform  flow  part 


?n  =  U x  -  r d  y . 


3.1.39) 


3.1.2  Inner  Solution 

The  inner  solution  will  be  obtained  by  stretching  the  coordinates 
with  a  suitable  stretching  factor,  in  this  case  the  reverse  of 
the  blade  aspect  ratio  t.  The  three-dimensional  Laplace  equa¬ 
tion  becomes  the  two-dimensional  one  with  an  accuracy  to  the 
second  order  of  .  The  boundary  conditions  near  and  on  the 
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blade  remain  the  same  whereas  those  at  infinity  points  lose 
their  meanings  due  to  the  stretching  of  the  coordinates.  Since 
much  literature  has  already  described  the  detailed  procedure 
in  similar  problems,  the  formal  coordinate  expansions,  etc., 
will  be  deleted  herein. 

Figure  3.1.2  shows  a  flow  configuration  for  the  inner  solution 
on  the  new  coordinate  system  attached  to  the  blade.  The  rela¬ 
tionship  between  the  old  and  new  coordinate  systems  is  given 
as  follows 


x  =  x  sin  B(r)  -  y  cos  6(r) 
y  =  -x  cos  ?  (r)  -  y  sin  3  (r)  . 


(3.1.40) 


or 

,  •  -*  i3 

Y  =  -i y  e 

where  _  _ 

Y  =  x  +  iy  ,  y  =  x  iy  ,  y  *  =  x  -  iy .  (3.1.41) 

The  selection  of  a  ca-  ity  closure  model  to  be  used  for  the 
present  singular  perturbation  problem  needs  particular  atten¬ 
tion.  In  this  problem  the  assumption  was  made  that  the  cavity 
length  be  short  and  can  collapse  to  a  line  in  the  coordinate 
stretching  procedure.  This  assumption  immediately  eliminates 
the  chance  of  using  a  cavity  closure  model  in  which  the 
information  in  the  inner  region  carried  over  to  the  outer 
region.  The  open  wake  model  and  double  spiral  vortex  model 
are  included  in  such  category.  The  only  feasible  one  is  the 
single  spiral  vortex  model  of  Tulin  (1964).  In  this  model 
the  physical  quantities  on  the  upper  and  lower  wakes  are 
assumed  to  be  identical  for  the  same  velocity  potential  value 
$  so  that  the  body-cavity  streamlines  are  closed  at  the  cavity 
end  point  in  the  complex  potential  plane. 

The  nonlinear  solution  method  for  supercavitating  flows  with 
the  single  spiral  vortex  model  has  already  been  described  in 
Section  2.2. 
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In  order  to  be  able  to  match  the  inner  solution  with  the  outer 
solution,  the  former  must  be  expanded  in  terms  of  the  velocity 
potential  t.  From  Equation  (2.2.3), 


dW  iw 
dy  =  —  e 


(3.1.42) 


but  dW/d'  can  be  calculated  from  Equation  (2.2.2),  thus 
iw 


dy  =  % 


2  2; 

3  '?  :  - ’ 


4  2,2 

a  (;  +  a  ) 


d; 


3.1.43) 


Expansion  will  be  carried  out  as  |  y  j -*■  »  or  ;-*-ia  in  the  '-plane, 
therefore . 

,  .  .2 

giw  (  5)  _  eijw(ia)+w,(ia)(?-ia)+co"(ia)  ■  --2— -a —  +  •  •  • } 


=  eljj(ia)  |l  +  K1-  (;  -  ia)  +  K2  -  (;  -  ia) 


■}  (3.1.44) 


where 


=  i  w 1 (ia) 


K2  = 


i  jj  1  ( ia) 
qi 

^(ia)  =  iT  +  i  in  — 

1  q„ 


-  |w’  (ia)}  2J  j 2. 


(3.1.45) 


Uj,  q^  =  incoming  flow  angle  and  velocity,  respectively,  at 
upstream  infinity  in  the  inner  flow  region. 


It  must  be  mentioned  that  i and  q  are  quantities  yet  to  be 
known  and  thus  should  not  be  interpreted  as  those  of  the  outer 
flow  region  at  this  stage  of  the  theory.  They  will  shortly  be 
determined  through  the  matching  procedure.  dy  in  Equation 
(3.1.43)  now  becomes 


qie 


-iaT 


ia : 


dy  ~ 


(  ;  -  i  a ) 


+  k2  +  0  - ia) 


as  C  ia  , 
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by  integration 


q^e 


lag 


-  _  -ia  +  Kl  l n(;  -  ia)  +  KQ  +  K 2(r  -  ia) 


+  0 


where  Kn  is  an  integration  constant. 


-  ia)  2J> 


(3.1 


Inversion  of  Equation  (3.1.46)  yields  to 


,  -  ,  -  -K  L  +K 

——7—  AnY  -K,  in  y  +  (-K,  in  A.+  K.)  +  K.  -- 

-  ia  0  1  1  0  0  1  A  t  Afty 


+  0  l^1 


~2  )  , 


where  Ag  =  q^-e  xa^/ia :  J2 


h 


L0  =  -Kx  inA0-K0 


which  is  substituted  into  W  in  Equation  (2.2.2),  providing 


-i^  -  ao ; 

W  ~  qIe  1  y  +  — y1  ^  ‘  (ia)  inY 


iaoi  ( 

+  — —  |  i 


(ia)  inAfl+K0j 


i  a<?  , 


+  |  w  ’  ( i  a ) }  — — - T 


qxe 


iai  Y 


.  -  2J,2  K.L  -K,  . 

in  v  +  a  l  _  1  0  2  1 


qie 


-iai  Y 


(■¥} 


as  Y 


(3.1 


Take  the  complex  conjugate  on  both  sides  of  Equation  (3.1.47) 
and  express  the  quantities  of  complex  conjugate  by  * , 


i  ot  T — * 

W*  —  q  e  iY  + 


a  5  , 


*  in  Y* 


~2~  +  M0  +  M1  T* 


+  M2*  4,  +  o  [%1 


(3.1 


i 


Y 


iap 

where  Mq  =  — 


l-  {-i,’ 


( ia)  in  Aq  -  Kg^ 


.46) 


.47) 


.48  ) 
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3.1.3 


Matching  Procedure 

The  matching  procedure  for  the  singular  perturbation  method 
is  also  explained  in  much  literature.  This  includes  the  works 
by  Van  Dyke  (1964),  Ogilvie  (1970),  Brockett  (1972)  and  James 
(1975).  It  is  stated  in  the  book  by  Van  Dyke  (1964)  that  the 
asymptotic  matching  principle  will  be  satisfied  by  "matching 
the  m-term  inner  expansion  of  the  n-term  outer  expansion  with 
the  n-term  outer  expansion  of  the  m-term  inner  expansion"  where 
m  and  n  are  any  two  integers.  Due  to  the  ample  literature  exist¬ 
ing  on  the  subject,  the  matching  steps  herein  will  be  followed 
without  formal  procedure.  It  was  for  the  reason  that  the 
inner  and  outer  expansions  of  the  outer  and  inner  solutions 
in  the  previous  sections  were  made  without  expanding  the  co¬ 
ordinate  variables.  However,  in  using  such  an  informal  approach, 
care  must  be  taken  for  the  order  of  magnitude  of  each  term  in 
matching.  The  matching  will  be  carried  out  for  the  velocity 
potentials  p  in  Equations  (3.1.38)  and  (3.1.52). 


The  zeroth  order  matching  will  determine  the  unknown  velocity 
amplitude  and  direction,  and  a  of  the  inner  flow  field  as 
follows : 


q^1)(r)  =  ju2  +  (rd  2j  2  +  0(e) 


(3.1.53) 


/  (1)  7T  x  riu  .  -1/t  „\ 

tan (&!  -t  +  T-)  =  —  =  tan  ^  ~ 


(r)  =  0  +  0(e) 


(3.1.54) 


where 


0(e)  <  0(1)  and 


AR  = 


Propeller  Blade  Area 
Propeller  Radius 
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It  should  be  mentioned  that  the  inner  flow  region  has  been 
selected  to  have  th  .  coordinate  system,  (x,y)  in  line  with  the 
hydrodynamic  flow  direction  made  with  U  and  wr  (see  Figure 
3.1.2).  The  superscripts  on  qx  and  oc  denote  the  quantities 
of  first  o:  der  matching. 

The  first  order  matching  will  further  be  carried  out  for  deter¬ 
mining  the  strength  of  circulation  As  a  result  of  the 

zeroth  order  matching,  the  inner  cavity  flow  problem  can  be 
solved  for  the  first  time.  The  boundary  conditions  available 
for  the  inner  problem  include  the  following: 


i)  At  upstream  infinity,  i.e.,  at  ;  =  ia, 

qU/ 

w(ia)  =  ot^^  +  iln  1  /qa 


(3.1.55) 


where  w  (; )  ,  a^,  and  q-j.  are  given  in  Equations  (2.2.4), 
(3.1.53)  and  (3.1.54),  respectively  and 


qc  =  qlU)  V1  +  a  ' 


ii) 

I 


The  closure  condition  (see  p.  173  of  the  paper  by 
Larock  and  Street  (1965))  is  given  by 


where  C  denotes  the  contour  enclosing  the  body-cavity 
system,  which  is  transformed  into  c-plane, 


(3.1.56) 
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with  the  residue  theorem  employed,  finally 

Re  joj’  (ia)J  =  0.  (3.1.57) 

where  w'  denotes  the  derivative  of  jj  with  respect 
to  £.  More  detailed  derivation  of  Equation  (3.1.55) 
from  Equation  (3.1.54)  is.  discussed  in  the  paper  by 
Furuya  and  Maekawa  (1930). 


iii)  The  scaling  condition  between  the  physical  and  trans¬ 
form  plane  states  that  the  wetted  portion  of  the  body 
calculated  from  the  theory  should  be  equal  to  S 
specified  in  the  physical  plane.  On  the  wetted  part 
of  the  streamline 


dY  - 

ds 


or 


ds 


(3.1.58) 


where  3 


tan 


-1 


(dy/dx) 


and  y  denotes  the  body  coordinate. 


For  -1<£  <b,  w(£)  can  be  written  as  follows: 


j  i  g(?)  +  3(1)  ,  -1  <£  <  0 

w(>)  =  { 

(  i  g  U)  +  8  (  £)  +TT,  0  <  £(  <  b 


where 

g(U=x/d+  C)(b(1)  -  0 
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waere 


indicates  a  Cauchy  principle  integral 
script  1  denotes  the  first  order  quantities, 
grating  Equation  (3.1.58)  we  obtain 


and  super- 
By  inte- 


s  ( ?) 


where 
sg  (?) 


?  >  0 

?  <  0. 


The  arc  length  condition  is  therefore  satisfied  by 
the  following  equation: 


S  -  s(-l)  =  0.  (3.1.59) 

The  four  unknown  parameters  in  the  inner  flow,  i.e.,  a(1  , 
bd)^  p(l)  an(j  ^  ^)  win  now  be  uniquely  determined  by  using 
the  above  four  boundary  conditions.  It  is,  of  course,  under¬ 
stood  that  these  four  equations  are  highly  retarded,  thus 
requiring  a  numerical  method  such  as  that  introduced  by 
Furuya  (1975)  for  a  similar  problem. 


Once  the  inner  problem  is  solved  above,  the  first  order  match¬ 
ing  for  F  is  readily  achieved  by  comparing  the  corresponding 
terms  in  f>  for  the  inner  and  outer  expansion  of  Equations 
(3.1.38)  and  (3.1.52),  i.e.. 


P(1)(r)  = 


a(D .  d) 
a  .p 


m  L 


r.'( 


ia 


(1) 


+  0(e) 


(3.1.60) 


It  must  be  mentioned  that  the  right-hand  side  of  Equation 
(3.1.60)  is  the  totally  known  quantity  obtained  by  solving 
the  inner  solution.  It  is  also  important  to  notice  for  the 
latter  matching  procedure  that  F(1)(r),  therefore  apn/2' 
Imfw'(ia)j  are  quantities  of  order  e. 
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The  source  term  SQ(r)  in  the  inner  expansion  of  the  outer  solu¬ 
tion  is  matched  with  the  corresponding  term  in  the  counterpart, 
but  it  should  be  zero  due  to  the  closure  condition  applied  to 
the  inner  solution  (i.e.,  Equation  (3.1.57)), 

(1)  (1) 

SQ(1)(r)  =  - | -  Re  [w'  (ia(1))J  =  0.  (3.1.61) 


In  order  to  carry  out  the  second  order  matching,  the  second 
order  inner  solution  must  be  considered.  As  has  been  mentioned 
in  the  beginning  of  Section  3.1.2,  the  two-dimensional  Laplace 
equation  still  holds  with  exactly  the  same  boundary  conditions 
and  thus  the  expansion  of  the  inner  solution  remains  the  same 
as  that  of  Equation  (3.1.51).  The  difference  in  matching, 
however,  exists,  particularly  in  the  outer  solution;  the  first 
order  matching  has  provided  the  value  for  7^  (r)  in  Equation 
(3.1.38).  Once  the  circulation  distribution  r^(r)  on  the 

lifting  line  is  known,  the  induced  velocities  u (1)(r)  and 

(1)  ^ 
ufc  (r)  in  Equation  (3.1.38)  will  be  readily  calculated  from 

the  formula  given  in  Equations  (3.1.17)  and  (3.1.19).  Vari¬ 
ous  methods  of  computing  these  integrals  are  available.  One 
of  the  most  popular  and  convenient  methods  will  be  that  by 
Lerb  (1952)  who  applied  Nicholson's  asymptotic  formula  for 
singular  integrals.  The  second  order  matching  will  first  be 
made  for  the  incoming  flow  velocity  and  direction; 


q 


(2) 

I 


U  +  u 


(1) 


(r) 


0(,2) 


(3.1.62) 


and 


m  rui  -  u  (1)  (r) 

tan  (a  -  3  +  T)  =  - - -  =  tan 

1  1  U  +  u  L  (r) 

a 


or 


i^2'>  =  +  0(e2) 


(3.1.63) 
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where 


a  ■ 

1 


tan 


-1 


rui  -  u 


U  +  u 


(1) 

t 

ITT 


a 


(r) 

(r) 


(3.1.64) 


and  the  induced  flow  angle  is  depicted  in  Figure  3.1.2. 


The  circulation  matched  to  the  second  order  will  now  be 
written  in  the  following  form, 

r(2)(r)  =  -  4  Im[w’(ia(2))]+  0(e3)  (3.1.65) 


where  all  quantities  of  the  right-hand  side  in  Equation  (3.1.65) 
must  be  obtained  with  a  new  set  of  the  upstream  flow  conditions, 
i.e.,  q^25  and  a  which  are  defined  in  Equations  (3.1.62) 
and  (3.1.63).  The  matching  for  the  source  term  remains  the 
same  as  before, 


.(2)  „  (2) 


Rep  (  ia(2))]  s  o  • 


(3.1.66) 


The  moment  terms  in  Equation  (3 
corresponding  terms  in  Equation 


.1.38)  will 
(3.1.50); 

0  (e3)  J 

O  U3)  \ 


be  matched  with  the 


(3.1.67) 


The  result  of  this  matching  indicates  that  the  moment  due  to 
the  circulation  distribution  or  that  due  to  the  source  distri¬ 
bution  cannot  be  distinguished  within  the  framework  of  the 
singularity  distribution  method.  However,  practically  more 
control  for  the  circulation  distribution  can  be  possible  than 
that  for  the  source  distribution  sc  that  the  former  will  be 
exercised  in  actual  propeller  design  work. 
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It  can  be  realized  now  that  the  results  of  the  matching  pro¬ 
cedure  take  a  somewhat  different  form  from  that  in  other  con¬ 
ventional  methods  such  as  those  of  Van  Dyke  (1964),  Ogilvie 
(1970),  Brockett  (1972)  and  James  (1975).  In  those  works, 
the  circulation  7  was  expanded  as  an  ascending  series  of  t, 


where 

This  type  of  expansion  was  not  possible  in  the  present  approach 
r.nce  the  nonlinear  theory  was  applied  to  the  inner  solution. 
Change  of  the  circulation  1?  due  to  that  of  the  upstream  flow 
condition  is  not  additive  to  7^.  Therefore,  after  each 
matching  procedure,  7  should  be  successively  rewritten  as 
follows 

7  (r)  =  7(1)(r)  +  O(s^)  after  the  first  order  matching 

(  2  )  3 

7(r)  =  7  (r)  +  0(s  )  after  the  second  order  matching 


and  so  on. 

Finally,  it  should  be  mentioned  that  the  third  term  from  the 

last  in  Equation  (3.1.50)  was  not  left  out  for  matching.  It 

is  readily  seen  from  Equations  (3.1.49)  and  (3.1.51)  that  M' 

is  of  an  order  e  so  that  the  whole  term  has  an  order  of 
2 

e  ?n  c.  This  is  a  higher  order  than  that  of  the  next  term 
and,  as  a  matter  of  fact,  the  same  order  as  that  of  the  last 
term.  Matching  for  this  term  cannot  be  possible  within  the 
second  order  matching  but  will  require  the  higher  order  inner 
and  outer  solutions. 
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3.1.4 


Propeller  Performance 

The  thrust  and  torque  coefficients  of  a  propeller  can  be  cal¬ 
culated  based  on  the  circulation,  flow  velocity  and  angle. 

The  local  lifting  force  AL  is  obtained  from  the  Kutta- 
Joukowski  law 


-Mr) 


U  +  u 


(1) 

a 


(r) 


(3.1.68) 


It  should  be  mentioned  that  all  quantities  in  this  equation 
are  dimensional  ones,  whereas  throughout  the  analysis  in  Section 
3  the  normalized  ones  were  used.  This  point  was  mentioned  in 
Section  2  in  which  the  nondimensional  f  was  defined  as 


u  (1)(r)  =  u  (1)(r)/U 

3.  cL 

ut(1)  (r)  =  ufc(1)(r)/U 
~(r)  =  T/2  UR 


r  =  r/R 


but  the  bars  above  the  letters  had  been  dropped  for  convenience. 
The  thrust  force  T^  and  power  P^  of  the  propeller  are  thus 
written 


Th  =  K  J'  AL(r)  cos(S(r)  +  ou  )  dr 


K 

P^  =  K  J'  ur  AL(r)  sin(8(r)  a/)  dr. 


Conventional  normalization  provides  the  thrust  and  power  coef¬ 
ficients,  C„  and  C  , 


* 


3.2  LONG  CAVITY  PROPELLER  CASE  (CASE  2) 

This  is  the  case  in  which  the  cavity  length  l  has  the  same 
order  of  magnitude  as  that  of  the  propeller  radius  R.  As  is 
shown  in  Table  3.0.1,  i.e., 

IT  -  0(1) 

whereas  the  chord  c  to  the  propeller  radius  ratio  remains  the 
same  as  before, 


c 

R 


0(e) 


AR 


_1_ 

AR 

Propeller  Blade  Area 
Propeller  Radius 


It  is  interpreted  that  when  the  chord  shrinks  to  a  lifting  line 
with  the  propeller  radius  (or  span)  fixed,  the  cavity  cannot 
shrink  to  a  line.  The  outer  solution  should  therefore  consist 
of  a  lifting  line  followed  by  a  cavity  sheet  as  is  shown  in 
Figure  3.2.1. (a) . 

The  inner  solution  obtained  with  x  =  :X,  y  =  eY  and  z  =  Z  as 
e  -*•  0  will  be  the  one  having  the  cyan  length  and  cavity  length 
to  be  out  of  sight.  Therefore,  the  formal  singular  perturbation 
method  would  request  the  inner  flow  configuration  to  have  an 
infinite  cavity.  However,  a  difficult  problem  of  using  such  a 
flow  configuration  arises  in  the  matching  procedure.  Since  the 
cavitation  number  a  may  be  defined  as 


where 


p  ,  p  =  static  pressures  inside  the  cavity  and  at  the 
upstream  infinity,  respectively, 
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=  velocity  at  the  upstream  infinity. 


the  velocity  on  the  cavity  wall  should  be  calculated  by  Bernoulli 
equation 


qc  V"'-  +  '  . 

For  the  infinite  cavity  problem,  the  upstream  flow  velocity  in 
the  inner  region  q^.  should  be  identical  to  q  .  As  a  result  it 
is  found  that 


qT  =  qc  =  q^V1  +  c-  (3.2.2) 

As  has  been  seen  from  the  previous  section,  the  first  order 
matching  will  provide 

qT  =  q„  0.2.3) 

which  conflicts  the  precedent  result  obtained  in  Equation  (3.2.3). 
This  problem  apparently  stemmed  from  the  erroneous  choice  of 
the  inner  flow  configuration.  The  rule  of  the  singular  perturba¬ 
tion  method  says  that  the  inner  or  outer  flow  solution  shall  not 
carry  any  physical  or  flow  information  connecting  the  counter¬ 
part  solution.  In  the  above  discussion  the  cavitation  number 
which  defines  the  physical  relationship  between  the  cavity  of 
the  inner  region  pressure  and  the  upstream-infinity  pressure 
at  the  outer  region  also  determined  the  upstream  pressure  or 
velocity  q  of  the  inner  region.  This  overspecification  for 
the  boundary  conditions  caused  a  problem  of  proper  matching 
for  determining  the  unknown  quantities. 

Because  of  the  reason  just  mentioned,  the  inner  solution  must 
be  obtained  by  solving  a  problem  for  the  flow  configuration 
having  a  finite  cavity  length.  Although  the  cavity  treated 
in  the  present  case  is  assumed  to  be  too  long  to  be  enclosed 
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in  the  inner  region,  the  complete  body-cavity  system  will  have 
to  be  squeezed  within  the  inner  region.  A  question  arises  as 
to  what  kind  of  closure  condition  is  to  be  used.  Due  to  over¬ 
stretching  the  inner  coordinate  in  the  present  case,  one  might 
consider  imposing  a  condition  that  the  total  source  of  the  body- 
cavity  system  be  finite  in  order  to  represent  the  long  cavity 
not  to  be  enclosed  in  this  region.  However,  there  is  no  way 
at  this  stage  to  determine  a  quantitative  number  for  the 
finite  total  source  strength.  It  is  for  this  reason  that  the 

first  order  inner  solution  will  use  S  =0  and  that,  if  there 

o 

is  any  correction,  it  should  come  from  the  matching  procedure. 


The  first-order  inner  solution  in  this  case  is  therefore  iden¬ 
tical  to  that  used  in  the  previous  section.  The  outer  expansion 
of  the  inner  solution  in  Equation  (3.1.52)  can  be  used  for  the 
matching  purpose. 


As  far  as  the  inner  expansion  of  the  outer  solution  is  concerned, 
the  lifting-line  circulation  and  its  free  vortex  have  the  same 
expansion  form  as  Case  1,  i.e.. 


v 


,  c ,  ( r )  •  x  +  t .  ( r )  ■  y 

(r)~  7(r)  tan  (u  -x+u.-y)  +  — -  +  H.O.T. 

X  d  C  / 


2  2 
x  +  v 


(3.2.4) 


The  difference  in  the  present  case,  however,  exists  for  the 
expansion  of  i  due  to  the  source  terms,  $  and  its  matching 
with  the  inner  solution.  As  has  been  mentioned  in  the  beginning 
of  this  section,  when  z  ■*  0  with  the  span  length  fixed,  the 
cavity  fails  to  shrink  to  a  line  in  the  present  case.  The 
cavity  sheet  is  left  over  behind  the  lifting  line.  Any  quan¬ 
tities  associated  with  this  cavity  sheet  cannot  be  determined 
by  matching  on  the  line  because  they  do  not  belong  to  the  same 
stretching  factor  e.  In  order  to  determine  s(r,0)  in  Equation 
(2.1.8),  the  matching  with  a  new  scaling  parameter  5  will  be 
necessary.  The  new  scaling  parameter  3  will  be  defined  as  the 


! 
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j 

ji 


:] 


maximum  cavity  thickness  to  the  cavity  length.  In  this  sense 
the  present  problem  is  categorized  as  a  "multiple"  scale  problem 
in  the  singular  perturbation  method  as  is  explained  in  the  book 
by  Van  Dyke  (1964). 


The  same  potential  function  due  to  the  source  singularity  as  before 
(Equation  (2.1.8))  is  used 


1  -E 


£  f  r  iitUili.  nut 

k=l  y  r  R~  • 


(3.2.5) 


rh  3, 


R 


a,  =  T(  x  -  A  3  ' 

=  b 


)  2  +  r  '  2  +  r  2  -  2r'r  cos  'f,  , 


AS  ' )  2  +  (r  ’  -  r  cos  4' a  , )  2  +  (r  sin  ’f  a  , )  2 


3.2.6) 


i  ,  ,  =  5  '  -  3  +  a, 
k. 


(3.2.  7) 


The  inner  expansion  of  the  outer  solution  for  the  new  parameter 
6  will  be  made  as  follows.  In  Equation  (3.2.5)  (s,r,9)  is 

expanded  as  a  point  of  (x,r,3)  approaches  the  blade-cavity  sur¬ 
face  in  the  direction  normal  to  its  surface  element.  By  changing 
the  coordinate  system  from  the  cylindrical  one  to  a  local  one 
(x,y,r),  attached  to  the  cavity-body  surface  as  is  shown  in 
Figure  3.2.2,  i>s  can  be  written 

r+Ar  x+Ax 


s5 


~*/  / 


s  (r ' , x  ’ )  dr  '  dx ' 


r-lr  5f-Ax 


(x-x')2  +  (r  -  r  ’)  *  +  y 


-21  4 


+H.O.T.  (3.2.3) 


i 


where 

s(r'  ,5  '  (x'))  =  s  ( r  '  ,  x  ' )  ^Jr  ' 2  +  •  c 
A  =  U/wR, 


* 


and  consideration  has  been  made  that  the  major  contribution  for 
p  -  arises  from  the  nearest  singularity  with  j.  =  0.  The  subscriot 
5  for  p  has  been  used  to  distinguish  the  expansion  in  terms  of 
5  from  that  in  ; . 


In  order  to  expand  p  »  as  y  ■+  0 ,  the  Fourier  transform  method 
used  in  Section  3.1.1  will  again  be  employed  (see  also  page  23 
of  Ogilvie  (1970)). 


Define  an  integral  p= ;  as  follows 

,x)  dr'  dx' 


>--\f  hr?. 7- 


—  CX>  —DO 


r-  -  2  2—27  — 

l(x  -  x')  +  (r  -  r')  +  y  2 


where 


s(r',x)  =  0  for  r'-r  _>  Ir 

,  x  '  -  x  ,  >  lx . 


(3.2.9) 


The  Fourier  transform  of  p  5  with  respect  to  x,  is  given  by 


1  f  C  f  s  ( r  '  ,  x  ' )  e  dr  '  dx  1  dx 

2J  J  J  n-  — , .  2  ,  ,.2,-214 

_°o  -oo  -oo  (x  -  x  )  +  (r-r')  +y  2 


With  a  change  of  variables,  t  =  x-x'  and  defining  the  Fourier 

transform  of  j  ,  and  s  by  o  .*  and  s*,  respectively, 

SO  so 


»  oo  DO 


s  : 


(k 


•  •  ',/If 


s  ( r  1  ,  x  1 )  e~xk(t  +  x']  dr,  dx,  dt 


—  OO  — oo  — oo 


|t2  +  (r  -  r  ' )  2  +  y2"1  " 


$ 


♦ 


OO 


oo 


s*  ( r  '  ,  k )  dr 


e 


-ikt 


dt 


(r-r') 


2 


1 

2 
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X. 


s* (r ' , k) dr  '  •  KQ 


2 


where  is  the  modified  Bessel  function  (see  page  9  of  Fourier 
Cosine  Transforms  Section  of  the  book  by  Erdelyi  (1954))  for 
the  half  integral  of  the  real  part.  The  second  Fourier  trans¬ 
form  is  taken  with  respect  to  r; 


J>  ( k  , m ; y ) 

S  j 


-If 

—  co  —co 


s* ( r  '  , k)  dr  1  e 


-imr 


CO  oo 

--If 


s* (r  '  ,  k)  dr  '  e 


,  -  im  ( t  +  r  ' )  T,  j 


Kft{  k  (t2  +  y2)2!dt 


—  oo  —oo 


=  -  S**  ( 


co 

(k,m)  J" 


-imt 


Kn{ ! k  !  (t2  +  y2) 2  ^  dt 


_  s *  *  ( k ,  m )  yj(k2+m2)2 

_  -tt  2  2  - 

k  +  m 


(3. 2 . 10  ) 


where  the  formula  for  the  Fourier  cosine  transform  of  KQ  on  page 
56  of  Erdelyi  (1954)  has  been  used  and  the  double  stars  **  for 
the  superscripts  denote  the  double  Fourier  transform  of  each 
function. 


Expansion  of  $s?**  for  small  y  can  now  be  made; 


+>s  5**  (k  ,  m ;  y)  =  -  t 


s**  (k,m) 

2  2  — 
(kz  +  m  ) 2 


1  -|y! (k2  +  m2)2  + 


1  -2 
—  x  r 

2  2 


2\ 

m  ) 


+  0 


(y3) 


]■ 


Inverting  Equation  (3.2.11)  term  by  term, 


(3.2.11) 


51 


5  ( x ,  r ,  y )  =  ?  (x,r,0)+n-  y  s  ( r ,  x ) 


1  -i  2 

+  2  17  sy ' 


[j>s5  (x,r,i 


0)  +  3  .  ( x , r ,  0 ) 

XX  S'  ' 


z]  *  0  (y: 


)  (3.2.12 ) 


where  $  , (x, r, 0)  and  ?  .(x,r,0)  denote  the  second  derivatives 
S  XX  so  z  z 

of  i  ,(x,0,z)  with  respect  to  x  and  z,  respectively  and 

1  2 

1  s  (r  '  ,  x ' )  dr  '  dx ' 

5  .  ( x ,  r ,  0  )=-•=•/  /  - 1 - p  (3.2.13) 

3  2y  7  r-  -n2. ,  m21t 

rh  o  jjx  -  x')  +  (r  -  r  ■  )  J2  . 


There  exists  a  problem  in  inverting  the  above  equation:  the 
first  term  in  Equation  (3.2.12)  f  (x,r,0)  cannot  be  expressed 

S  1  j 

directly  in  terms  of  s(r,x)  due  to  the  extra  term  in  Equation 

2  2  —  —  — 
(3.2.11),  i.e.,  (k  +m  )2.  In  order  to  determine  s(r,x),  the 

second  term  must  be  used  as  will  be  seen  shortly. 

The  inner  solution  to  be  matched  with  the  outer  solution  expan¬ 
sion  in  Equation  (3.2.13)  should  be  obtained  from  stretching 
the  coordinates  based  on  the  two  parameters,  e  and  6.  With 
x  =  £ X ,  y  =  £ Y ,  z  =  Z  and  5  =  maximum  cavity  thickness/  ic 
as  £  and  :  -0,  the  flow  configuration  will  become  the  one  used 
in  a  linearized  supercavitating  flow  theory;  the  body-cavity 
streamlines  are  mapped  onto  a  thin  slit  as  is  shown  in  Figure 
3.2.1(b).  Although  it  may  be  possible  to  expand  the  nonlinear 
solution  around  the  y-axis,  such  expansion  seems  extremely  com¬ 
plicated  due  to  the  nature  of  the  theory.  The  purpose  of  using 
the  linearized  theory  here  is  simply  to  determine  the  source 
strength  for  the  cavity  thickness.  The  determination  of  all 
other  major  characteristics  still  depends  on  the  nonlinear 
theory.  The  higher  accuracy  of  the  overall  results  will  thus 
be  maintained  with  the  major  portion  including  the  circulation 
to  be  determined  by  the  nonlinear  theory. 

As  has  been  described  in  Section  2.3  for  the  2-D  linearized 
supercavitating  flow  theory,  the  source  distribution  on  the 
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body-cavity  slit  can  be  obtained  as  a  function  of  the  cavita¬ 
tion  number  and  geometric  foil  profile.  The  velocity  potential 
due  to  the  source  distribution,  expanded  in  terms  of  small  y, 
will  be  expressed  in  the  following  form 


* / 

S 


V 


where 


v 


(x) 

2 

(x) 

2 


0  <  x  <-1 

1  <-  x  <2 

c 


(3.2.14) 


(3.2.15) 


and  s^(x)  and  s2(x)  have  been  given  in  Equations  (2.3.16)  and 
(2.3.17). 


Matching  between  the  inner  and  outer  solutions  will  now  deter¬ 
mine  the  unknown  source  distribution  function  s(r,x); 


(x) 

2-r 

s2  (x) 

s ^ ( x )  in  Equation  (2.3.16), 
s2 (x)  in  Equation  (2.3.17). 


s (r, x)  = 


0  <  x  <1 

1  <  x  <  l 

c 


(3.2.16) 


Once  the  source  distribution  is  obtained,  the  outer  solution 
due  to  the  souce  singularity  will  be  expanded  now  in  terms  of  e. 
As  the  blade  chord  shrinks  to  a  line,  is  expressed  in  the 
similar  form  to  that  in  Equation  (3.1.11), 
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K 


k=l 


x  sn  (r  1  )  +  r  sin?.  •  s,  (r  '  ) 
lx _  k  ly _ 

R„ 


dr 


I  AI'  I 


s_  (r  1  ,6  ■  (x1)) 


d5  '  +  H.O.T. 


(3.2.17) 


where 


SQ(r)  =  J  s1(r',e’  (x’))  d3  ' 


(3.2.18) 


Slx  {r 


')  =  J  s1(r',S'  (x* ))  S'  (-Xr* 


)d9  ’ 


(3.2.19) 


-f 


s.  (r ' )  =  /  s.  (r ’ ,9 ' (x')  0  '  d6 ' 


j 

'L 


(3.2.20) 


and  and  s2  are  given  in  Equations  (2.3.16)  and  (2.3.17). 
The  last  term  in  Equation  (3.2.17)  is  the  new  term  evolving 
from  the  existence  of  the  long  cavity  sheet. 


Now  the  inner  expansion  of  the  outer  solution  as  c  -*•  0  becomes 


~Sn(r)  ?.n  (x  +  y 


sc  0 


2  .  ..2,f  Sln(r)  'x  *  Slv(r>  ,y 


2  .  2 
x  +  y 


+  (“as’ 


+  H.O.T. 


(3.2.21) 


where 


3 


s2 


as 


3  x 
3  J 


u 


s2 


ts  r3  9 


(3.2.22) 
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(3.2.23) 


The  two-term  inner  expansion  of  the  outer  solution  j  is  now 
written 

(ua  ( r )  +  uas  (r))-x 


*>0  +  (r)  tan 


_1 *  +  S  (r) in (x2  +  y2) 2  + 
x  0  1 


Uix(r)-x+  Uly-y 
2  X  2 

x  +  y 


+  H.O.T. 


(3.2.24) 


where 


$  =  Ux  -  rcu  y . 


(3.2.25) 


The  same  matching  procedure  as  in  Case  I  can  be  applied  here. 

This  will  result  in  almost  the  identical  solution  as  before 

except  for  the  incoming  flow  velocity  qT,  induced  velocities 

( 2 )  1 

ou,  and  the  source  strength  ; 


(2) 


♦  r  -  u'1’ <r,  -  u<l>  <r> 


1 

21  2 


+  0(e‘ 


(3.2.26) 


(2)  2 
=  a  +  0(e 


I 


ot .  =  tan 
l 


-1 


ru)  -  u 


U  +  u 


(1) 

t 

nr 


(r)  -  ufcs  (r) 


(r)  +  u  (r) 

a,  S 


-(H- 


(3.2.27) 

(3.2.28) 


(2) 


(2)  (2) 


a  p 


i 


Re  '  (ia  ^  2  ^  )J  =  Sq  in  Equation  (3.2.18). 


(3.2.29) 
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The  result  of  the  second  order  matching  indicates  that  the 

second  order  inner  problem  to  be  solved  has  the  incoming  flow 

(2)  (2)  ( 2 )  — 

quantities  q^  and  with  the  closure  condition  Sq  =  SQ, 

finite.  It  means  that  (1)  the  upstream  flow  now  has  the 

induced  velocity  correction  not  only  due  to  the  vortex  wake 

but  also  the  source  sheet  and  (2)  the  total  source  strength 

is  not  equal  to  zero,  but  finite,  given  by  Equation  (3.2.18). 

The  latter  condition  is  interpreted  as  a  correction  to  the 

assumption  made  for  the  first  order  inner  solution,  as  has 

been  mentioned  earlier.  The  matching  for  the  circulation  7 

of  the  second  order  provides  exactly  the  same  form  as  that  in 

Case  I ,  i . e . , 


J2)(2) 

’2)  ?l 

'  J  (r)  =  - 5-^-  I 

2  m 


[<*>’  (iaC2))]  _ 


(3.2.30) 


Table  3.2.1  is  provided  to  help  clarify  the  present  matching 
procedure  in  comparison  with  Case  I. 


4.0  CONCLUSION  AND  RECOMMENDATION 

Unlike  finite  span  wings  or  subcavitating  propellers,  the 
singular  perturbation  method  for  the  supercavitating  propeller 
required  careful  classification  of  the  problem  due  to  the 
existence  of  the  multiple  scaling  parameters.  The  scaling 
parameters  here  included  1)  span  length  R,  2)  chord  length 
c,  3)  blade  spacing  d  and  4)  cavity  length  1  .  The  span 
length  R  was  chosen  to  be  the  reference  parameter  for  con¬ 
structing  various  scaling  parameters,  namely  c/R,  d/R  and  /.  /R. 

The  first  problem  solved  here  assumed  that  c/R  and  l  /R 
were  of  order  e  but  d/R  was  of  order  of  unity.  It  turned 
out  that  the  nature  of  the  singular  perturbation  problem 
was  similar  to  that  for  the  subcavitating  propeller  solved 
by  Brockett  (1972)  except  for  the  solution  of  the  inner 
region.  The  thrust  and  torque  coefficients  were  obtained 
explicitly  without  solving  the  integral  equations.  Since  the 
nonlinear  supercavitating  flow  theory  was  employed  in  the 
present  work  as  the  inner  solution,  there  existed  no  limitation 
for  the  flow  incidence  angles  or  blade  profile  shapes.  Due 
to  the  nature  of  the  nonlinear  theory,  the  loading  coefficients 
could  not  be  expressed  in  ascending  series  of  t,  but  were  cal¬ 
culated  with  new  boundary  conditions  applied.  It  must  be  men¬ 
tioned  that  the  new  calculation  will  not  cause  any  difficulty 
since  the  same  formula  for  the  boundary  value  problem  can  be 
utilized  with  changing  only  the  boundary  conditions.  It  is 
believed  that  the  present  solution  obtained  with  the  nonlinear 
theory  as  the  inner  solution  will  provide  more  accurate  results 
than  those  with  the  linearized  theory. 

The  second  problem  treated  in  this  study  was  the  case  in  which 
c/R  had  an  order  of  e  but  £c/R  and  d/R  were  of  order  of  unity. 

Physically,  this  was  the  case  having  long  cavities  behind  the 
propeller  blades  so  that  even  when  the  chord  shrank  to  a  line, 
the  cavities  were  left  behind  the  lifting  lines.  This  portion 

l 
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of  cavity  sheets  was  called  "source  sheets",  the  singularity 
strengths  of  which  were  obtained  through  the  cavity  sheet 
matching.  This  matching  was  totally  different  from  the  regular 
matching  carried  out  for  the  properties  belonging  to  the  lift¬ 
ing  line.  It  was  considered  that  the  present  problem  had  to 
be  categorized  as  the  multiple  scale  problem  in  the  singular 
perturbation  method. 

The  first-order  inner  solution  used  a  closure  condition, 
i.e.,  the  total  source  term  equal  to  zero.  It  was  con¬ 
sidered  that  this  assumption  might  not  be  totally  correct, 
at  least  physically,  because  the  cavity  length  was  too  long 
to  be  fully  contained  in  the  inner  region.  As  the  result  of 
the  cavity  sheet  matching,  however,  it  was  discovered  that 
the  second-order  inner  solution  had  to  use  a  finite  SQ,  the 
value  of  which  was  determined  through  the  matching  procedure. 

It  seemed  that  the  second-order  matching  automatically  cor¬ 
rected  the  overstretching  assumption  made  in  the  first-order 
inner  solution. 

In  the  section  of  problem  classification,  Section  3.0,  two 
other  problems  were  posed,  both  having  small  blade  spacings. 

The  outer  solutions  may  be  quite  different  from  those  in 
the  above  two  cases  since  if  one  looks  at  such  a  propeller 
from  the  far  field,  the  blade  elements  will  not  be  identified. 
This  may  require  the  actuator  disc  concept  for  the  outer  solu¬ 
tion,  with  a  conventional  pressure  jump  across  the  disc  if  the 
cavity  is  short  but  with  a  cavity  pressure  drop  to  be  applied 
if  the  cavity  is  long,  similar  to  the  theory  of  Tulin  (1965). 
These  two  problems  have  not  been  carried  out  here  due  to  the 
enormous  amount  of  work  required  even  for  the  first  two  problems. 

For  the  future  study,  therefore,  the  last  two  supercavitating 
problems  having  a  large  number  of  blades  are  recommended  to 
be  solved  with  the  singular  perturbation  method.  It  is  also 
our  regret  that  numerical  computations  have  not  been  conducted 


58 


for  comparison  of  the  present  analytical  results  with  experi¬ 
mental  data  for  the  same  reason  above.  Such  comparison  will 
be  interesting  and  is  also  recommended  for  further  efforts. 
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CLASSIFICATION  OF  S/C  PROPELLER  PROBLEMS  IN  SRM 


TABLE  3.2.1 


MATCHING  PROCEDURE  AND  RESULTS  FOR  CASE  I  AND  CASE  II 


CASE  I  (Short  Cavity) 


1st  Order 

2nd  Order 

a{0) 

. (0)  _ ( 1 ) 

/ 

S'**' 

Outer 

Solution 

V 

,  (2) 

-  (2) 

3o  = 

Inner 

Solution 

(1) 

q 

(1)  „  (2)  (2)  ^ 
r  1  'if1 

-  (1) 

=  0  s = o 

CASE  II  (Long  Cavity. 

.  '  ‘..s-  :  r  :  2nd  Order 


Outer  Solution 


Inner  Solution 


(a)  Physical  plane  7  =  x  +  iy 


in 

k 

ia 


b  d 


(c)  Transform  plane  ^  =  -;  +  in 

Flow  configuration  and  transform  planes 
for  the  same  supercavitating  flow  as  that 
in  Ficure  2.2.1  but  with  double  spiral 
vortex  model 


FIGURE  2.2.2 


.ate  system  (x,y) 
.he  blade  and 
■  conditions 


APPENDIX  A 

INVERSION  FORMULA  OF  SINGULAR  INTEGRAL  EQUATIONS 


The  inversion  formula  of  the  singular  integral  equation  for 
arcs  is  well  described  in  the  book  of  Muskhelishvili  (1946). 

In  order  to  provide  a  brief  insight  into  the  derivation  of  the 
formula  the  inversion  method  will  be  summarized  herein.  Let 
the  singular  integral  equation  be  defined  by 


1  f  ?  ( t)  dt 

t~t0 

L 


f ( tg )  on  L 


( A-l ) 


where  the  contour  L  consists  oi 


smooth  arcs  L,  to  L  , 
1  P 


i .  e .  , 


L  —  L^  +  L2  +  *  •  •  +  L^,  (A— 2) 

f(t)  is  a  given  function  and  i(t)  is  to  be  determined.  The 

functions  f(t)  and  d(t)  are  considered  to  belong  to  the  classes 
1)2) 

H  and  H*  ,  respectively. 


Note:  1)  A  function  ?(t)  will  be  said  to  satisfy  a  Holder  condi¬ 

tion  (or  H  condition)  on  the  arc  L,  if  for  any  two  points 
tx,  t2  of  L 


|5>(t. 


-  <j)  (t1)  I  1 A  s  t2  -  t1 


where  A  and  u  are  positive  constants.  The  function  ?(t) 
will  be  said  to  belong  to  the  class  H  on  L,  if  it  satisfies 
the  H(h)  condition  for  some  u  >  0  on  each  of  the  closed 
arcs  Lj  of  L  including  the  ends. 


2)  If  the  function  $(t),  given  on  L,  satisfies  the  H(y) 
condition  on  every  closed  part  of  L  not  containing  ends, 
and  if  near  any  end  c  it  is  of  the  form 


1  \ ot  = 

( t  -  c) 

where  i*(t)  belongs  to  the  class  H,  then  t(t)  will  be 
said  to  belong  to  the  class  H*  on  L. 


73 


Introducing  a  sectionallv  holomorphic  function 


MZ) 


1  f  D  (t)dt 
2tt1  J  t  -  Z 
L 


(A-3) 


one  has 


* 

V 


+ 


(t. 


)  +  5 


(t0) 


1  f  2  (  t)  dt 

Ti  J  t-  tQ  • 

L 


Therefore,  the  inversion  of  the  singular  integral  equation  (A-l) 
is  equivalent  to  the  problem 


$+(tQ)  +  ( tp )  =  f(tQ)  on  L,  and 

$(»)  =  0. 


(A-4) 


This  is  a  well-known  mixed-type  boundary  value  problem,  the 
solution  method  for  which  requires  first  the  homogeneous  solu¬ 
tion  X(Z)  of  ( tg )  +'S”(tCJ)  =  0.  By  choosing  a  function 


X  (Z) 


R2 ( Z )  >2 


(A-5) 


where 


r,  (z)  =  n  (z  -  c.) 

1  k=l  K 

2p 

r,(z)  =  n  (z  -  c,  ) 

k=q+l  K 


( A-6 ) 


( A-7 ) 


and  the  quantity  of  X(Z)  is  understood  to  refer  to  the  branch 
cut  along  L.  It  is  readily  found  then  that 


[x  ( Z  )J  =  -  jx  ( z  )j 


(A-8 ) 


and  therefore  a  solution  to  the  problem  for  a  new  function  Y(Z] 
Y  ( Z  )  =  ^  (Z)/x  (Z)  (A-9  ) 
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is  given  by  solving  a  boundary  value  problem 


Y 


+ 


-  Y 


t>(z)]+  +  r$(z)i 

[X(Z)] + 


f  (t) 
[X(t)]  + 


on  L 


( A-10 ) 


The  solution  for  Y(Z)  is  given 


Y  (Z) 


1  f  f  (t) 

2Ti  J  [X(t)]  + 

L 


Q  (Z) 


(A-ll ) 


where  Q(Z)  is  an  arbitrary  polyncminal.  The  solution  for  ;(Z) 
is  therefore  obtained  by  substituting  (A-ll)  into  (A-9) , 


3  (Z) 


X(Z) 


f  (t) 
[X(t)]  + 


dt 
i  -  Z 


+ 


Q(Z) 


(A- 12 ) 


The  arbitrary  polynominal  Q(Z)  will  be  chosen  in  such  a  way 
that  J>(=°)  =  0  and  other  physical  conditions  of  the  problem 
such  as  singularities  be  satisfied. 
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